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We discuss a multichannel SU(JV) Kondo model which displays non-trivial zero-temperature phase 
transitions due to a conduction electron density of states vanishing with a power law at the Fermi 
level. In a particular large-TV limit, the system is described by coupled integral equations corre- 
sponding to a dynamic saddle point. We exactly determine the universal low-energy behavior of 
spectral densities at the scale-invariant fixed points, obtain anomalous exponents, and compute scal- 
ing functions describing the crossover near the quantum-critical points. We argue that our findings 
are relevant to recent experiments on impurity-doped d-wave superconductors. 



The Kondo effect in metals is a well studied phe- 
nomenon in many-body physics. The low-energy physics 
is completely determined by a single energy scale, the 
Kondo temperature 7k, and the impurity spin is fully 
quenched in the low-temperature limit, T <C 7k The 
standard picture of the Kondo effect has to be revised 
if the conduction band density of states (DOS) vanishes 
identically at the Fermi energy This is the case 

in so-called pseudogap systems with a power-law DOS 
p(e) ~ |e| r (r > 0) which arises in one-dimensional in- 
teracting systems, in certain zero-gap semiconductors, 
and in systems with long-range order where the order 
parameter has nodes at the Fermi surface, e.g., p- and d- 
wave superconductors (r = 2 and 1). Of special interest 
are high-T c cuprate superconductors, where indeed non- 
trivial Kondo-like behavior has been observed associated 
with the magnetic moments induced by non-magnetic Zn 
impurities 0. 

The pseudogap Kondo problem has attracted a lot of 
attention during the last decade. A number of studies 
||,|| including the initial work by Withoff and Frad- 
kin employed a slave-boson l&rge-N technique; further 
progress and insight came from numerical renormaliza- 
tion group (NRG) calculations ^-^] and the local mo- 
ment approach ^ . The general picture arising from these 
studies is that there exists a zero-temperature phase 
transition at a critical Kondo coupling, J c , below which 
the impurity spin is unscreened even at lowest temper- 
atures. Also, the behavior depends sensitively on the 
presence or absence of particle-hole (p-h) symmetry: in 
the p-h symmetric case there is no complete screening 
even for 7k > 7 C . A comprehensive discussion of possi- 
ble fixed points and their thermodynamic properties has 
been given by Gonzalez-Buxton and Ingersent |(| based 
on the NRG approach. Spectral properties of the pseu- 
dogap Kondo model have so far been investigated for the 
p-h symmetric case only H[|, an d the knowledge about 
the quantum-critical dynamics is limited. Importantly, 
recent advances in scanning tunneling microscopy (STM) 
have made it possible to directly observe local spectral 
properties of impurities in superconductors || . 

The purpose of this paper is to provide an exact low- 



energy solution of the pseudogap Kondo problem in a cer- 
tain large- TV limit [|l0|]ll| l which is appropriate for obtain- 
ing dynamical properties. We shall determine anoma- 
lous exponents and spectral densities at the fixed points 
of interest, and describe the crossovers in the vicinity 
of quantum-critical points. The critical fixed points are 
found to have exponents varying continuously with the 
pseudogap exponent r. To our knowledge, our study rep- 
resents the first critical theory of the pseudogap Kondo 
transition in the presence of p-h asymmetry. 

We begin by describing the generalized multichannel 
Kondo problem, characterized by K channels of conduc- 
tion electrons and a spin symmetry group extended from 
SU(2) to SU(7V). The choice of the SU(TV) representa- 
tion for the impurity spin S qualitatively influences the 
phase diagram in the large- TV limit. For r = 0, a symmet- 
ric (bosonic) representation leads to both underscreened 
and overscreened |l^,13| phases with an exactly screened 
phase in between pC|, whereas an antisymmetric repre- 
sentation only allows for an overscreened situation f pT| . 
A second crucial difference is that the bosonic large-TV 
limit contains a generic p-h asymmetry which prevents a 
study of a p-h symmetric model. Therefore, we choose an 
antisymmetric representation [ [PH for S, corresponding to 
a Young tableau with a single column of Q boxes. The 
TV 2 — 1 traceless components of S can be written in terms 
of TV auxiliary fermions, /„, as S ufl = flf^ — QS^^/N, 
together with the constraint Ylufvfv — Q- The con- 
duction electrons transform under the fundamental rep- 
resentation of SU(TV), and carry a channel (flavor) index 
i = 1, • • • , K and SU(TV) spin index v = 1, •••,7V. The 
Hamiltonian of the model reads: 
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The fermionic bath is assumed to have a power law DOS 
at low energies; for simplicity we use a model DOS given 
by p(e) = Pa \e/D\ r for |e| < D. 

For r = 0, it is known that weak antiferromagnetic 
coupling 7k grows under renormalization for K > 2 and 
any TV, Q, and there is a stable intermediate-coupling 
non-Fermi-liquid fixed point associated with overscreen- 
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FIG. 1. Schematic phase diagrams of the multichannel 
pseudogap Kondo model as function of the Kondo coupling 
Jk and the particle- hole asymmetry (e.g. potential scatter- 
ing), for pseudogap exponents r with a) < r < r* lin , b) 
r mm < r < r* mx for 7 > 7 C (top) and 7 < y c (bottom), 
7c ~ 0.8, c) r > r* lax . 

ing of the impurity spin. For r > 0, we find that the 
situation is different, and a number of phases appear: (i) 
A local moment (LM) regime is reached for weak initial 
coupling. Here Jk flows to zero, and the impurity is un- 
screened, (ii) An overscreened phase (OS) is found for 
< r < f*Jn ax and J c i < Jk < J C 2 , it is a natural general- 
ization of the r = overscreened fixed point, (iii) There 
is a new asymmetric strong-coupling (ASC) phase which 
exists for r > 0, and is reached only for Jk > J C 2 and a 
sufficient amount of p-h symmetry breaking. Qualitative 
phase diagrams are shown Fig. |l|; the important case of 
r = 1 corresponds to Fig. ^3. 

We now proceed with the analysis of the model in the 
large- iV limit which has been introduced in Refs. p0| , pl| . 
We set K = 7 TV, J K = J/N, and take the limit N -> 00 
keeping 7 and J fixed. The important difference of the 
present large- N limit to the usual non-crossing approx- 
imation (NCA) is the treatment of the occupation 
constraint: here we set Q = qoN and keep qo fixed. The 
advantage is twofold: we can achieve p-h symmetry with 
Q = N/2, and in the large- N limit the dynamics of the 
system is determined by a true (time-dependent) saddle 
point ]l(j[n|. (Formally, results corresponding to Q = 1 
as in Rcf. |14[] are obtained here in the limit qa — > 0.) The 
Kondo interaction in Eq. ([!]) can be decoupled in each 
channel by auxiliary bosonic fields Bi(r), conjugate to 
the amplitude ^ u c\ v {t)J v (t). The constraint is imple- 
mented with a Lagrange multiplier field A(r). One finds 
the following integral equations for the fermionic and 
bosonic Green's functions G/(r) = — (T r / v (r)/ J J(0)), 
G b {t) = <T t £? 4 (t)/4(0)): 

S / (r)= 7 G (r)c7 B (r) , E s (r) = -G Q {-r)G f (r) , (2) 
with the self-energies E f and E b defined by: 

Gj 1 (iu} n ) = iuj n + A - Hf(iu n ) (3) 



G^ilVn) = — - E.b(m/„ 



J 



(4) 



The third saddle point equation fixes the impurity "oc- 
cupation" and thus determines A: 

G/(r = 0-) = ^G f (iu> n )e^ 0+ = q . (5) 

n 

Eqs. (j|,||,||) are identical in structure to the usual NCA 
equations, but the physics is changed by the constraint 
cq. (||) which keeps track of the choice of the impurity 
spin representation. With a symmetric conduction band 
DOS, the model is p-h symmetric (invariant under /J «-> 
f u and c\_ w <-> c kw ) for q = 1/2. P-h symmetry can 
be broken at the impurity spin (qo ^ 1/2) as well as in 
the environment (potential scattering or doping), both 
leading to similar effects. 

An important dynamic observable is the conduction 
electron T-matrix T(cu) (which corresponds to the impu- 
rity Green's function for an impurity Anderson model), 
given by the convolution of G/ and Gb, so its Fourier 
transform is T(r) = -(l/JV)G/(r)G B (-r). Local sus- 
ceptibilities for the spin and channel degrees of free- 
dom are obtained by xi oc ,sp = — Jq drG / (r)G / (— t) and 
Xioc,ch = Jq rfTGs(T)Gs(— t); these have to be distin- 
guished from the impurity contributions Ximp,sp,ch to the 
reponsc to a global (spin or channel) field. 

We now turn to a detailed analysis of the low-energy, 
low-temperature behavior of the above equations. In 
a first step, we analyze possibly existing scale-invariant 
fixed points where the Green's functions follow power 
laws at small frequencies. It is important to note that 
any continuous quantum phase transition (as well as 
other possible intermediate coupling fixed points) must 
be described by such behavior. The constant terms 
in Eqs. (§,|) have to vanish, A - E/(i0 + ) -» 0,1/ J - 
E B (i0 + ) — > as T — > fn|]. A power-law ansatz 
G b \t) - A b /t 2Ab , G/(r) ~ A f /r 2A f for the Green's 
functions in the regime 1/T*<§;t</3^oo (where T* 
is a crossover temperature) is equivalent to low-energy 
spectral densities 

( 6 ) 



which have the same power- law exponent, but different 
prefactors A±, for ui > and uj < 0. Inserting these 
Green's functions into Eqs. (|^J|,(§) one finds the following 
equations for the exponents: 



2A/ + 2AB = l-r, 
(1 - 2A s )sin27rAij _ (1 - 2A/) sin27rA / 



sin ttAb + sinh a/2 cos 2 it Af + sinh a/2 



(7) 



Here, a is a spectral asymmetry parameter defined by 



exp(a) = Af t —/Af t . 



-A B -/A B .+ . For r = 0, a 



In these expressions uj n = (2n + 1)tt//3 and v n = 2m:/ (3 
denote fermionic and bosonic Matsubara frequencies. 



drops out, and there is a single solution with non-trivial 
exponents, 2A/ = 1/(1 + 7), 2Ab = 7/(1 + 7), which 
describes an infrared stable fixed point corresponding to 
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FIG. 2. Values of the exponent 2A/ at the scale- invariant 
fixed points, obtained from Eq. ([?]), as function of r. Solid: in- 
frared stable fixed points (OS). Dashed: unstable fixed points 
(quantum phase transitions). Left: particle-hole symmetric 
case a = 0. Right: with particle-hole asymmetry. The disap- 
pearance of the stable OS solution marks the value of r*(a). 
We find r min , r^ ax -> 1 for 7 -> 0, r min ~ 0.19, r max ~ 0.29 for 
7=1, and r mm , r maJ[ — * for 7 — > 00. (In the single-channel 
SU(2) case, r* = §, whereas N = K = 2 gives r* ~ 0.23 ||.) 

an overscreened spin pl] |. In contrast, the behavior for 
< r < 1 is much richer, and is illustrated in Figs. |l| and 

For small r and a^fl, one finds three solutions, one 
infrared stable representing the OS phase, and two unsta- 
ble corresponding to the quantum phase transitions from 
OS to LM and ASC. With increasing r, two solutions dis- 
appear at r* (a) , leaving a single unstable solution for the 
direct transition between the LM and ASC phases. In the 
p-h symmetric case (a=0), the unstable solution for the 
transition to ASC is absent. Therefore, no scale-invarian 
solution is found for r > r*(0): the p-h symmetric sys- 
tem here always flows to the LM weak coupling phase. 
The critical value r*(a) is bounded, r^ in < r* < r^ ax . 
From the presence or absence of intermediate coupling 
fixed points one can infer phase diagrams as in Fig. [l|; 
these are also well borne out by our numerical studies of 
Eqs. (|!||i|). I n the narrow range r* lin < r < r^ ax , the 
existence of the OS phase depends on the amount of p-h 
symmetry breaking, Fig. ^b. 

The N = K = 2 version of the overscreened pseudogap 
Kondo model has been studied earlier by NRG ||, with 
results somewhat different from the large- TV limit. The 
N = K = 2 model has a single r* value, no ASC region for 
r < r* ~ 0.23, and a minimal potential scattering value 
to reach ASC for r > r* . Some of these differences arise 
from the fact that the N = K = 2 case exhibits isolated 
critical fixed points, whereas the scale-invariant large- N 
solutions actually correspond to continuous lines of fixed 
points, parametrized by the low-energy spectral asym- 
metry a, which is determined by the symmetry-breaking 
parameters of the model. 

The universal low-energy behavior of the T-matrix 
at the scale-invariant fixed points is readily obtained 
from the convolution of G/ and Gb , its spectral density 
Pt(uj — > 0) diverges as AtM -5 ". Notably, the prefactors 
for both signs of uj are equal here, i.e., the asymmetry 



described by a disappears in the scaling limit Fur- 
thermore, the amplitude At in the OS phase is indepen- 
dent of J, this is the analog of Friedel's sum rule. The 
scaling forms of the finite-temperature spectral densities 
can be determined [jll] from the T = quantities (||) by 
applying a conformal transformation z = exp(i27rr//3). 
Low-temperature thermodynamic quantities can now be 
calculated, e.g., the local spin susceptibility (l^Jl^] fol- 
lows xioc, S p ~ l/T 1 "^/ for 4A/ < 1, X ioc,s P ~ ln(l/T) 
for 4Ay = 1, and xi oc ,sp ~ const otherwise; similar re- 
sults hold for xioc,ch in terms of Af,. 

In a second step, we discuss possible non-scale- 
invariant fixed points of the system (which are absent 
at r = 0): the local moment (LM) and the asymmetric 
strong coupling (ASC) fixed point. These solutions of 
Eqs. (HJl,i|) appear if the T = "compensation" of the 
constant terms in Eqs. (||^) is not realized. 

For the LM fixed point, reached for small Kondo cou- 
pling, it is important to notice that Es(iQ+) is bounded 
from above for r > and cannot compensate large 1/J 
values (the corresponding integral is infrared divergent 
for r = 0). Hence, there is no power-law solution for 
J smaller than a critical J ci . A constant term in G^ 1 
implies a <5-peak in the T = fermion spectral density; 
furthermore the spectral asymmetry flows to zero in the 
scaling limit, leading to pf(uj) ~ 5(u>). This is equivalent 
to an unscreened spin: the local spin susceptibility fol- 
lows a Curie law for T — > 0. The "local moment" Txioc.sp 
is reduced from its free spin value, but Txi mPjSp , measur- 
ing the total impurity-induced Curie contribution to the 
uniform susceptibility, is that of a free spin, i.e., part 
of the moment is carried by the surrounding conduction 
electrons due to the residual coupling mediated by Gb- 

The ASC fixed point is reached at sufficient particle- 
hole asymmetry and Kondo coupling. Here, A — S/(i0 + ) 
does not vanish for T — > 0. This leads to a <5-peak in the 
auxiliary boson spectral density, ps(w) ~ ±5(w— 0^) [the 
sign depending on sgn(go — 1/2) ], which implies maximal 
spectral asymmetry in the scaling limit. The channel 
susceptibility diverges like l/T, and the system has a free 
"moment" in the channel (flavor) degrees of freedom, also 
leading to a residual T = entropy. In contrast, the spin 
is completely screened, Xim P , sp , Xioc, sp — > const. The 
ASC phase can be understood as (N — Q) conduction 
electrons bound to the impurity; for N = K = 2 this is 
equivalent to a spin singlet, channel doublet state ||. 

For both the LM and ASC fixed points, the spectral 
density pt{uj) of the T-matrix behaves as \u\ r at small 
frequencies, and the spectral asymmetry disappears in 
the scaling limit of px- 

At this point we briefly comment on the case r > 1: 
the OS phase disappears for all 7; in the presence of 
particle-hole asymmetry a transition between LM and 
ASC is still possible. For r = 1 the above power law 
analysis remains valid, but has to be supplemented with 
logarithmic corrections. For r > 1 the critical power law 
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FIG. 3. Zero-temperature scaling function $t for the con- 
duction electron T-matrix at r — 0.8, 7 = 1, qo = 0.1, describ- 
ing the crossover from the ASC behavior at low energies to 
quantum criticality at high energies. Inset: Same on log-log 
scale, showing both u> > (solid) and ui < (dashed). 

at the LM-ASC transition is replaced by /?t ~ M r ~ 2 - 

After having described the fixed points we now turn 
to the universal behavior near the quantum phase tran- 
sitions. In the vicinity of each critical point, one can de- 
fine an energy scale T* , which vanishes at the transition, 
and defines the crossover energy above which quantum- 
critical behavior is observed |15| . From the above low- 
energy analysis one can deduce the dependence of T* on 
the reduced coupling, j = (J — J c )/J c , measuring the 
distance from the critical point at J c . Approaching the 
LM-OS or LM-ASC transition from small J we have T* ~ 
{-j) 1/r , from large J wc find T* ~ f/(r+2& f ,cr) wne re 
Ay cr is the anomalous fermion exponent at the critical 
fixed point. Similarly, we have T* ~ (_j)i/(r-+2A / , os ) 



and T* 



,-l/(r+2A /iCr ) 



if we approach the OS-ASC 



transition from below or above, respectively. The lo- 
cal susceptibilities are easily derived, e.g., in the ASC 
phase we have xioc. sp (T=0) ~ T*( 4A /,cr-i) _ j-i w j t h 
7= (l-4A / , cr )/(r'+2A / , cr ). 

Near a critical point all observables show scaling be- 
havior as function of T/T* and to/T* (k B = 1). For 
instance, the T-matrix obeys the scaling form 

Tip) = r . (1 _„ T) *t f ^ , 7^ J , (8) 



where A is an amplitude prefactor, 7/t = 1 — r is the 
anomalous exponent, and <I>t is a universal scaling func- 
tion (for the particular transition and values of r, 7, 
go)- As an example, we focus on the transition be- 
tween the LM and ASC phases, for J > J c . A result 
for the T = scaling function $t is shown in Fig. [| For 
u> = uj/T* «lwe have ASC behavior with Im<I>T ~ |w| r , 
for w > 1 the spectral density follows the quantum- 
critical power law |ct>|~ r as predicted above. In both lim- 
its the spectrum is particle-hole symmetric (!), however, 
in the crossover region the asymmetry leads to a large 
peak for one sign of Q [depending again on sgn(<7 — 1/2) ]. 

Finally, we mention possible applications of our re- 
sults. Although we have taken a multichannel large-Af 
limit here, the T-matrix spectral properties of the LM, 



ASC, and critical fixed points turn out to apply also to 
the single-channel case [H. In most systems there is a 
finite particle-hole asymmetry, and if Kondo screening is 
present in a pseudogap system (and r is not too small), 
the physics will be described by the ASC fixed point M . 
Then, in contrast to the metallic r = case, the impu- 
rity spectral function does not show a peak at the Fermi 
level, but at an energy of order Tk ~ T* (Fig. [| see 
also ||). Interestingly, STM studies of Zn impurities 
in the cuprate <i-wave superconductor BSCCO || have 
shown a large peak in the differential conductance at a 
small bias of 1.5 meV, and a likely possibility is that this 
peak arises from the Kondo screening of the impurity- 
induced moment fl7j . The value of T* = 15 K is con- 
sistent with NMR studies which show a Tk of order 
20-40 K around optimal doping. Last not least, genuine 
multichannel pseudogap physics may be found in quan- 
tum dots coupled to d- wave-superconducting leads, and 
possibly in heavy- fermion superconductors [ fl3|| . 
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